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Seismic modeling of the Slowly Pulsating B-type star HD 182255 yields strong constraints on the radial orders of the two 
dominant modes. For these frequencies, we derive also the empirical values of the complex nonadiabatic parameter /. 
These two seismic tools, i.e., pulsational frequencies and the associated values of /, allow to test available opacity data, 
chemical composition and overshooting efficiency. 
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1 Introduction 

Slowly Pulsating B-type stars (SPB) pulsate in high order g- 
modes, which penetrate the deep interior of the star. These 
modes could be useful in particular to determine the over- 
shooting efficiency from the convective core. Because of a 
very dense oscillation spectrum of high order g-modes, seis- 
mic modeling the SPB stars is quite challenging if no reg- 
ular patterns are found. However, it appeared that fitting at 
least two frequencies reduces greatly possible combinations 
of the radial orders. Our analysis of the SPB star HD74560 
(Walczak, Szewczuk & Daszynska-Daszkiewicz 2012), has 
shown that only two combinations of the radial orders are 
possible. 

In our seismic modelling we include also the empiri- 
cal values of the nonadiabatic complex /-parameter, de- 
fined as the ratio of the radiative flux change to the radial 
displacement at the level of the photosphere (Daszynska- 
Daszkiewicz, Dziembowski & Pamyatnykh 2003, 2005). 
The values of / depend on the properties of the superficial 
layers and, in case of B-type stars, are sensitive to stellar 
opacities and metallicity. Fitting simultaneously frequen- 
cies and the /-parameters (so called complex asteroseismol- 
ogy) can provide better constraints on the parameters of the 
model and theory. 

In the next section, we give basic information about the 
star HD 182255. Results of the mode degree identification 
are presented in Section 13.11 In Section 13.21 we determine 
the radial orders of dominant modes and section |4]contains 
description of complex seismic models of the star. Conclu- 
sions end the paper. 
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2 HD182255 

HD 182255 is a bright star in the constellation of Vulpec- 
ula with about 5mag in the Johnson V filter. The spectral 
type of the star is B6 and the Hipparcos (ESA 1997) par- 
allax amounts to about 8.30 mas. According to Glebocki & 
Gnacinski (2005) the projection of its rotational velocity is 
equal to V smi =30 km/s. Hube & Aikman (1991) discov- 
ered that the star is in a binary system with the orbital pe- 
riod of about 367.26 days. A slightly larger period (367.76 
days) was found by Dukes et al. (2003). The later authors 
discovered also three other repeatable variations with peri- 
ods of ~ 1 day. Similar values have been derived by Math- 
ias et al. (2001). These variations were interpreted as pul- 
sations in high order g-modes. Five pulsational frequencies 
were found by De Cat et al. (2007) in data gained during 
photometric and spectroscopic campaign devoted to SPB 
stars (Aerts et al. 1999). All frequencies were present in the 
Geneva photometry while the two dominant ones were de- 
tected also in spectroscopy (De Cat et al. 2007). 

In order to estimate the basic parameters of the star, in 
Fig. Q] we show its position in the HR diagram. The ef- 
fective temperature, logT e ff = 4.149 ± 0.012, was taken 
from Niemczura (2003) who determined this value from the 
IUE (International Ultraviolet Explorer, Bogges et al. 1978) 
spectra. The luminosity, log L/Lq — 2.457 ± 0.087, was 
derived from the Hipparcos parallax and bolometric correc- 
tion from Code at al. (1976). We took into account also the 
interstellar extinction in the form of Ay = RE(B — V), 
where we adopted the standard value of R = 3.1 and the 
reddening of E(B - V) = 0.002 ± 0.011 (Niemczura 
2003). The evolutionary tracks from the Zero Age Main 
Sequence (ZAMS) to the Terminal Age Main Sequence 
(TAMS) were calculated by using the Warsaw-New Jer- 
sey evolutionary code. We adopted the AGSS09 chemical 
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mixture (Asplund et al. 2009) and the OPAL opacity tables 
(Iglesias & Rogers 1996) supplemented by the Ferguson, 
Alexander & Allard (2005) opacities. Furthermore, we as- 
sumed the initial rotational velocity of V TO t = 30 km/s, ini- 
tial hydrogen abundance X = 0.7, metallicity parameter 
Z = 0.02 and no overshooting from the convective core. 
The lines marked as g. L will be described in the next Sec- 
tion. The mass of HD 182255 estimated from evolutionary 
tracks is about 4M and most probably the star is in the 
main sequence phase. 

3 Mode identification 

3.1 The mode degree, t 

To determine the mode degree of the five detected fre- 
quencies of HD 182255 (De Cat et al. 2007), we compared 
firstly the theoretical values of the photometric amplitude 
ratios and phase differences with the observational coun- 
terparts using the theoretical values of the nonadiabatic /- 
parameter. The theoretical photometric observables were 
computed for models located in the center and on the edges 
of the observational error box. In the next step, the mode 
degrees were derived together with the empirical values of 
the /-parameter. This method was applied only to the two 
dominant frequencies because in case of the B-type stars 
the empirical value of the /-parameter can only be de- 
rived for frequencies which are visible in both photometry 
and spectroscopy (Daszyhska-Daszkiewicz, Dziembowski 
& Pamyatnykh 2005). The Geneva photometric system con- 
sists of seven filters and the response functions of the adja- 
cent ones overlap. Therefore, we used only UBV which can 
be treated as independent. The above two approaches need 
input from models of stellar atmospheres. We used the LTE 
models by Kurucz (2004) and the non-linear limb darkening 
law by Claret (2003). 

In Table [T] we present a summary of our identification 
of £. In the first two columns we give frequencies and the 
amplitudes of the light variation in the Geneva V filter, re- 
spectively. The most probable mode degrees derived with 
the first and second methods are given in the third and fourth 
columns, respectively. The first method relied on photome- 
try while the second - on photometry and spectroscopy. 

As we can see, from the first and second method we 
obtained that the dominant frequencies, v\ and are 
the dipole modes. The second method indicates also the 
quadrupole modes. The dipole modes are also the and 
v$ frequencies. The only ambiguous result is for the 1^3 fre- 
quency, which can be I — 1 or 2. 

Unfortunately, we do not have information about az- 
imuthal orders, m. Identification of m requires the analysis 
of line profile variations or taking into account rotational ef- 
fects which are beyond the scope of this paper. In the further 
analysis, we assume that pulsational modes of HD 182255 
are axisymmetric, i.e. m = 0. 



Table 1 Pulsational frequencies of HD 182255, the 
Geneva V amplitudes and the most probable identification 
of £ derived with the theoretical and empirical values of the 
/-parameter. 



Frequency 


A v 


Mode degree, £ 


[c/d] 


[mmag] 


theoretical / 


empirical / 


^1=0.97185(3) 


12.4(1.1) 


I = 1 


£ = 1,2 


^2=0.79225(4) 


13.4(1.0) 


I = 1 


£ = 1,2 


!/ 3 =0.62526(7) 


8.3(1.0) 


£=1,2 




^4=1.12780(10) 


3.5(1.0) 


I = 1 




^5=1.02884(11) 


2.4(1.1) 


I = 1 





3.2 The radial order, n 

We started our modeling by fitting the dominant frequency 
v\ = 0.97185 c/d assuming different radial orders, n. Fig.Q] 
shows the position of models fitting v\ in the HR diagram in 
cases when n has a value from 10 to 15. Each line is marked 
with the corresponding value of n. 
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Fig. 1 The position of HD 1 82255 in the HR diagram. Different 
symbols indicate models fitting the v\ = 0.97185 c/d frequency 
assuming radial orders from n = 10 to n = 15. All models and 
evolutionary tracks were calculated with the OPAL opacities and 
metallicity Z = 0.02. We did not include the overshooting effect. 



As we can see there are a lot of models at different evo- 
lutionary stages as well as various radial orders allowed. In 
general, for a given radial order we have a more evolved 
star when the stellar mass increases. And we get more mas- 
sive stars with decreasing radial order of the v\ mode. These 
effects are easily explained by the facts that for high-order 
g-modes the following rules apply: 

- for a given mass the frequency decreases with increasing 
radial order, 

- for a given radial order, the frequency increases with de- 
creasing stellar mass, 
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- for a given radial order and mass, the frequency gener- 
ally increases with time. 
When we try to fit a constant value of a frequency assuming 
increasing radial order we have to reduce the mass. With 
smaller masses we get lower effective temperatures and lu- 
minosities And for a given radial order, the bigger the mass 
we choose, the more evolved a star becomes. However, we 
have to bear in mind, that there are exceptions to these rules 
connected with the avoided crossing which appears for non- 
radial modes and modifies the frequency evolution. 

In the next step, we plotted these models in the fre- 
quency vs. mass diagrams. We chose the modes with fre- 
quencies closest to the observational value of V2 = 0.79225 
c/d. An example of such diagram in shown in Fig. [2] where 
we plotted the frequency of the £ = 1, g 15 mode as a func- 
tion of the mass. All these models fit the v\ frequency un- 
der assumption that it is the 1=1, gi2 mode. The hor- 
izontal line indicates a value of the v% frequency. As we 
can see, the frequency dependence of the mass is very com- 
plicated, but there are models fitting simultaneously the v\ 
and V2 frequencies. These models have the masses of about 
4.15M , 4.23M , 4.27M Q , 4.31M , 4.375M , 4.42M Q 
and 4.45M . 
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Fig. 2 Frequency of the 1=1, gis mode as a function of the 
mass for models fitting v\ (assuming that v\ is the 1 = 1, gi2 
mode). The results are for metallicity of Z = 0.02. A horizontal 
line indicates the value of the U2 frequency. 

In Fig. Q] we can easily see, that for Z = 0.02 some 
models calculated with the radial order n = 12 fit quite 
well the observational effective temperature and luminosity. 
We searched for models fitting two frequencies assuming 
different radial orders of the v\ frequency. As a result we 
found that models fitting v\ and v% and laying close to the 
observational error box (which means closer than the 3cr 
errors) exists only if v\ is the 1=1, gi2 mode and v-i is the 
1 = 1, gi5 mode. 

4 Complex asteroseismic analysis 

In this section we consider only seismic models fitting si- 
multaneously the two dominant frequencies: v\ {£ — 1, gi2) 



and v-i (£ = 1, gi5). Such models are plotted in the HR dia- 
gram in Fig. [3] where we show the metallicity effect. Along 
the lines connecting points there should also exist seismic 
models. However, between the lines there are no seismic 
models. The grey area indicates the models with the un- 
stable £ = 1, gi2 mode while the hatched regions indicate 
models with the unstable £ = 1, gis mode. The area labeled 
as f{v\) will be discussed later. 
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Fig. 3 The HR diagram with the location of the seismic mod- 
els fitting the two dominant frequencies: v\ [£ = 1, gi2) and vi 
{£ = 1, gis). Symbols mark the values of metallicity, Z. Gray and 
hatched regions indicate domains where v\ and U2 are unstable, 
respectively. The area labeled as f(ui) is described in the text. 



The effective temperature of seismic models is highly 
sensitive to the metallicity parameter. In general, a smaller 
value of Z gives higher effective temperature. We can also 
notice that for a given metallicity we have rather a large 
spread in luminosity while the effective temperature is only 
slightly variable. In order to catch models inside of the er- 
ror box we need quite a high metallicity parameter; from 
about Z = 0.02 up to 0.025. This is rather inconsistent with 
the metallicity derived by Niemczura (2003): [m/H] = 
-0.30 ±0.31 which corresponds to Z = 0.0098±[^Jg . All 
models lying inside of the error box have unstable modes 
corresponding to v\ and v-i . 

Because the two dominant frequencies were detected in 
both photometry and spectroscopy we were able to derive 
the empirical values of the non-adiabatic /-parameter and 
compare them with theoretical counterparts. Unfortunately, 
we did not fully succeed. We managed to find models fitting 
the /-parameter only for the v\ mode and within the 3a er- 
rors. The models lie inside of the area labeled as f(vi) in 
the HR diagram (Fig. [3J. It can be easily seen, that they are 
far from the error box; they have too high effective temper- 
ature. The /-parameter indicates very young models, close 
to ZAMS, and metallicity lower than 0.015. 
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One possible explanation of this disagreement between 
theoretical and empirical /-parameters is connected with 
the data. The problem is that, the spectroscopic observa- 
tions were made about two years before photometry. The 
value of the empirical /-parameter is sensitive to both pho- 
tometric and spectroscopic amplitudes and phases. If the 
amplitudes or phases changed significantly between obser- 
vations, the obtained value of / would be incorrect. This is 
justified suspicious since some seasonal variations in ampli- 
tudes and phases in Stromgren filters were found by Duces 
et al. (2003). 

We checked also whether our models fit other frequen- 
cies. It turned out that there is a model which reproduces 
quite accurately the 2/3 and v± frequencies. Its parameters 
are: M = 4.239M Q , logT cff = 4.1593, log L/L Q = 
2.4826 and Z = 0.02. In Fig. g] we plotted the instability 
parameter, 77, of this model as a function of frequency for 
modes with I = 1 — 4. Modes with 77 > are exited in the 
model. The vertical lines indicate observational frequencies 
of HD182255. The v z = 0.62526 c/d frequency can be both 
the dipole gig or quadrupole g 3 4 mode, but if it is the I — 2 
mode, it is stable (77 0.2). The 1/4 = 1.12780 c/d fre- 
quency is well reproduced by the dipole mode gio while 
the v$ = 1.02884 c/d frequency does not have a theoretical 
counterpart. The rotational splitting of the £ = 1, gis mode 
is about 0.2 c/d (assuming rotational velocity of 30 km/s) so 
7/5 and v\ could belong to the same triplet. It is interesting, 
that observational frequencies are inside of the instability 
domain of the dipole mode. 
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Fig. 4 Instability parameter as a function of frequency for 
modes with t = 1 — 4. The chosen model is described in the text. 
The vertical lines indicate observational frequencies. 



We took into account also overshooting from the con- 
vective core assuming a ov — 0.1, 0.2 and 0.4. Seismic mod- 
els fitting the v\ and v 2 frequencies changed slightly their 
position in the HR diagram. Moreover, we were not able 
to find models fitting other frequencies. Also the theoretical 
values of the non-adiabatic /-parameter did not fit better the 
empirical counterparts than in case of a ov = 0. 

Quite similar models appeared with the OP opacity ta- 
bles (Seaton 2005). The main difference was connected with 



the instability regions: with OP they were slightly larger. 
The quadrupole g34 mode nearly fitting the v z frequency 
was almost unstable for Z > 0.02. But still we could not 
find models fitting the empirical values of the /-parameter 
with effective temperatures and luminosities in agrement 
with the observational counterparts. 

5 Summary 

Fitting two observational frequencies allowed us to con- 
strain radial orders of the pulsational modes. Seismic mod- 
els close to the observational error box exist only for one 
combination of the radial orders: v\ is the I = 1, gi2 mode 
and V2 is the 1 = 1, gi5 mode. Including overshooting from 
the convective core did not improve our models. Moreover, 
with high values of a ov we could not find models fitting 
more than two frequencies. It seems, that effective over- 
shooting for HD 182255 is not necessary. 

We did not succeed in fitting the empirical and the- 
oretical values of the non-adiabatic /-parameter for both 
modes v\ and v-i. The possible explanation is the inaccu- 
racy in opacities on which the /-parameters strongly de- 
pend and/or lack of the simultaneous photometric and spec- 
troscopic observations. Further observations are needed as 
well as identification of the azimuthal number m for all pul- 
sational modes. 
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